Robust control of quantum gates via sequential convex programming 
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Resource tradeoffs can often be established by solving an appropriate robust optimization problem for a vari- 
ety of scenarios involving constraints on optimization variables and uncertainties. Using an approach based on 
sequential convex programming, we demonstrate that a substantial fidelity robustness is obtainable against un- 
certainties while simultaneously using limited resources of control amplitude and bandwidth. What is required 
is a specific knowledge of the range and character of the uncertainties, a process referred to in the control theory 
literature as "uncertainty modeling." Using a general one-qubit model for illustrative simulations of a controlled 
qubit, we generate robust controls for a universal gate set. Our results demonstrate that, even for this simple 
model, there exist a rich variety of control design possibilities. 



I. INTRODUCTION 

Robust control and robust optimization of uncertain sys- 
tems are essential in many areas of science and engineering 
[1-7], Recently, there has been much interest in achieving 
robust control of quantum information systems in the pres- 
ence of uncertainty [8-39]. An important property of quantum 
information processing that distinguishes it from most other 
applications is the requirement of an unprecedented degree 
of precision in controlling the system dynamics. Also, due 
to the very fast timescale of physical processes in the quan- 
tum realm, implementing closed-loop feedback control is ex- 
tremely difficult and thus open-loop control arises as the most 
feasible option in most circumstances. 

For quantum information systems, a robust optimization 
problem can be formulated as a search for design variables 
9 that maximize a measure of quantum gate fidelity T over a 
range of uncertain parameters 5 (i.e., (5 £ A, where A is the 
uncertainty set). Fidelity compares a target unitary transfor- 
mation against the actual quantum channel which is dependent 
on both 6 and S. Fidelity is typically normalized: < < 1, 
and the maximum value = 1 coiTesponds to a perfect gen- 
eration of the target transformation. The design variables 9 
can include time-dependent control fields (for both open-loop 
and closed-loop control), measurement configurations (usu- 
ally, for closed-loop feedback control), constants associated 
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with physical implementation, the circuit lay-out, and so on. 
The uncertainties in 5 can affect any elements of the system 
Hamiltonian (including the design variables), with specific 
manifestations and ranges depending on details of the phys- 
ical implementation and external hardware. For example, un- 
certainties can represent dispersion and/or slow time variation 
of parameters such as coupling strengths, exchange interac- 
tions, and applied electromagnetic fields, as well as additive 
and/or multiplicative noise in control fields. The uncertainty 
set A can thus, in general, combine deterministic and random 
variables. Whatever the case, we assume that both 9 and S are 
constrained to known sets 9 and A, respectively. 

One common approach to robust control of quantum gates 
(e.g., see Ref. [ ]) is based on maximizing the average fi- 
delity given by 

where J- {9, 5) denotes the fidelity as a function of design 
and uncertain variables, and E^gAl } is expectation with re- 
spect to the underlying distiibution in A. Often the aver- 
age fidelity is well approximated as the sum over a discrete 
sample and associated probabilities {Si e A,pi,Vi}, i.e., 
J'a.vg{9) = '}2nPiJ-{9 ,5i). While this approach is applicable 
in some cases (in particular, when the uncertainty represents 
weak random noise), the stringent performance requirements 
of quantum information processing make it more appropriate, 
in general, to estimate gate errors by using a worst-case fi- 
delity measure with respect to all uncertainties 5 G A, 

F^,{9)=minF{9,5). (2) 



Also, worst-case robust optimization (or minimax optimiza- 
tion) is a well known approach employed in many classical 
problems [7, 40-55], and some of the methods developed for 
these applications can be adapted for robust control of quan- 
tum gates. The worst-case robust optimization problem for 
quantum gate fidelity is formulated as: 

maximize xmnFiO^S) 

^ (3) 
subject to e 9, (5 e A 

The goal reflected in (3) is to find the design variables 6* e 
that maximize the worst-case fidelity (2). 

In control applications, the design set <d represents the set 
of control constraints. Q is most often convex or sufficiently 
well approximated by a convex set. In some cases so is the 
uncertainty set A, although even if it is not that is not the 
principal obstacle to solving (3). What makes the problem 
difficult is that the fidelity measure is not a convex function 
of 6 for any sample (5 € A. Non-convex optimization prob- 
lems are common in all of science and engineering and have 
engendered numerous numerical approaches to finding local 
optimal solutions. In particular, effective methods have been 
developed in recent years for worst-case robust optimization 
with non-convex cost functions [50-56]. 

In optimal control applications, the functional dependence 
of the objective (e.g., fidelity) on the control variables is re- 
ferred to as the optimal control landscape [57-60]. For an 
ideal model of a closed quantum system with no uncertain- 
ties, the optimal control landscape for the generation of uni- 
tary transformations has a very favorable topology. Specifi- 
cally, provided a number of physically reasonable conditions 
are satisfied, the landscape is free of local optima, i.e., there 
exist one manifold of global minimum solutions (resulting in 
F — Q) and one manifold of global maximum solutions (re- 
sulting in = 1), while all other extrema reside on saddle- 
point manifolds [61-63]. Such a favorable landscape topology 
facilitates easy optimization, as any gradient-based search is 
guaranteed to reach the global maximum [ ]. Unfortunately, 
when uncertainties are present, this landscape topology is not 
preserved. Typically, uncertainties cause the drop and frag- 
mentation of the global maximum manifold, resulting in the 
emergence of multiple local maxima [ ] (the landscape also 
undergoes a similar transformation when control fields are 
severely constrained [ ]). Provided that the range of uncer- 
tainty is not too large, many of these local optimal solutions 
will have fidelities close to one. 

For quantum information systems, there is considerable on- 
going effort to develop efficient methods for obtaining a good 



solution to the problem of robust control, for either average 
or worst-case fidelity. The majority of existing approaches 
rely on a numerical optimization procedure, mostly involving 
a gradient-based search for maximizing the average fidelity of 
Eq. (1). In some cases, a randomized search such as a genetic 
algorithm is employed [39]. The results demonstrate the ex- 
istence of many solutions with high fidelities, consistent with 
the control landscape picture discussed above. Additionally, 
the optimal controls are often similar to the corresponding ini- 
tial controls, provided the latter are reasonably good. This 
phenomenon, also observed in many engineering and design 
applications employing local search algorithms, supports the 
need for developing tools to efficiently calculate a good ini- 
tial control. In particular, empirical evidence and simulations 
suggest that robust controls for an uncertain quantum system 
can be found by searches that start from solutions generated 
by applying optimal control theory or dynamical decoupling 
to the ideal (zero-uncertainty) counterpart system (see, e.g., 
[60], [38] and the references therein). 

In this paper, we propose the use of sequential convex pro- 
gramming (SCP) which is one of several methods available 
for numerically solving optimization problems like (3). (See 
[ ] for a collection of earlier SCP varieties and uses and [67] 
for a recent informative overview.) SCP provides a frame- 
work for finding local optimal solutions to the general robust 
optimization problem (3). The specific SCP algorithm used 
here, delineated in (4) below, follows directly from [ , ]. 
It was used previously for robust design of slow-light tapers 
in photonic-crystal waveguides [54, 55] and quantum poten- 
tial profiles for electron transmission in semiconductor nan- 
odevices [56]. In this paper, we apply this SCP algorithm to 
identify robust control fields for generation of quantum gates 
in an uncertain one-qubit system. 

II. SEQUENTIAL CONVEX PROGRAMMING 

The SCP algorithm used here is shown in abstract form in 
(4). The algorithm must be initialized with (i) a control in the 
feasible set 8, which is assumed to be convex, (ii) samples 
5i,\ — 1, . . . ,L taken from the uncertainty set A, which need 
not be convex, and (iii) selecting the initial size of a convex 
trust region Btmst- The trust region is selected so that the lin- 
earized fidelity used in the optimization step retains sufficient 
accuracy. In each iteration the SCP algorithm returns the op- 
timal increment 6 and the associated worst-case (linearized) 
fidelity. To compute the actual worst-case fidelity requires 
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simulating the system with the control variables + 9 as in- 
dicated in Step 3 of (4). The centerpiece is the optimization 
step which, in the version shown (4) is gradient based, thereby 
resulting in L affine constraints in 9, and hence is a convex 
optimization. The Hessian, perhaps not so easily computed, is 
easily added as shown in Appendix A. In some cases the num- 
ber of samples, L, can be very large. Fortunately, however, 
computational complexity grows gracefully with the number 
of constraints and thus does not grossly effect the convex op- 
timization efficiency [68]. 



Robust control via SCP (4) 

Initialize 

initialize control 6* G 9 C 

sample uncertainties/noises 5i £ A, i = 1, . . . , L 
set size of trust region Strust CI R^ 

Repeat 

1. Calculate fidelities and gradients 
F{e, S,), VeHd, S,) eR"", i = l,...,L 

2. Using the linearized fidelity, solve fi>r the incre- 
ment 6 from the convex optimization: 

maximize min T{9, Si) + \/eJ~{9, Si)"^ 6 

i— 1,...,L 

subject to 9 + e e Q, 6 e Q trust 

3. Update 

If min T{e + e,5i)> min T{e,5i) 

i—l,...,L i — l....,L 

Then 9^9 + 6 

increase size of trust region Gtiust 
Else decrease size of trust region Otrust 

Until Stopping criteria satisfied 



In §A we show how the gradient and/or Hessian forms can 
be cast in standard forms compatible with freely available 
software specifically designed to solve such convex optimiza- 
tion problems. In general, solving the convex optimization 
step is not the most time consuming step in the SCP algorithm. 
The time-burden in each iteration falls more often on the sim- 
ulations required to obtain the fidelities and gradients (and the 
Hessian if used) at each uncertainty sample. Of course as is 
the case with numerically computing anything to do with a 
quantum information system, there always lurks the "curse of 
dimensionality", i.e., exponential scaling with the number of 
qubits. 



Despite many advantages, SCP is a local optimization 
method. As such there is no way to verify that an optimal 
solution has been found. In the case of maximizing fidelity, 
which by construction cannot exceed one, it would seem that 
at least the maximum is known, so if = 1 is achieved, it is 
an optimal solution. However, in many cases we see fidelities 
thatareextremely close to one, e.g., log ]^q(1 — J^) e [—6,-4]. 
Although 4 to 6 nines following the decimal point is 1 for 
most engineering problems, for quantum computing every 
additional nine can greatly decrease the spatial and tempo- 
ral resources required for fault-tolerant operation. Ergo, here 
Ml! 

III. SEQUENTIAL CONVEX PROGRAMMING FOR AN 
UNCERTAIN QUBIT 

In this section we show how to use SCP (4) for some com- 
mon uncertainty and control generation and constraint scenar- 
ios. We focus on controlling a single qubit system with Hamil- 
tonian. 

Hit) = c{t)uj^X + LJzZ, <t <T (5) 

where c{t) is the external control field and X, Z are the re- 
spective Pauli matrices. The parameters uJx and uj^ are con- 
stant but uncertain over the time interval [0, T]. In regard to 
(3), we can equate the uncertain parameters 5 with (w^;, Wz). 

A. Control generation and constraints 

The control c{t) is typically the output of a signal gener- 
ation device whose dynamics impose constraints on magni- 
tudes, bandwidth, and so on. To illustrate the use of SCP 
we make the simplifying assumption that the control is piece- 
wise-constant over N uniform time intervals of width T/N, 
i.e., for fc = 1, . . . , A^, 

c{t, 9) = 9k and (fc - l)T/N < t < kT/N, (6) 

Constraints due to signal generation dynamics are explored in 
§B. In general the set 6 in (3) reflects the control constraints. 
Consider, for example. Table I, showing many typical restric- 
tions on the control. 

This list of control constraints is certainly not exhaustive. 
However, since c(t, 9) is a linear function of 9 (6), then these, 
or any combinations, form convex constraints on9 £ . For 
example, with the strictly piece-wise-constant control (6) the 
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Constraint 


Set e 


none 




fluence 


/o^c2(t,e)dt<7 


magnitude 


c""° < c(t,e) < c'""", t£ [0,T] 


slew rate 


\c{t,e)\<l3, t€[0,T] 


area 


J^\c{t,e)\dt<a 


linear 


Ae = b 



TABLE I: Typical control constraints. 



fluence constraint becomes (r/A^)||6'||2 < 7^. The bounding 
parameters in Table I can also be used as design variables to 
establish control resource tradeoffs via SCR This is explored 
further in § V for fidelity vs. fluence. 



B. Fidelity 



For the single-qubit system (5), we use the fidelity, 



(7) 



This fidelity, normalized to [0, 1], is a measure of the phase- 
independent alignment between the desired unitary logic gate 
V and the actual unitary operator Ut attained at the final time 
T. Applying the piece-wise-constant control (6) to the Hamil- 
tonian (5) gives. 



the covariance matrix associated with a specified confidence 
region for the parameters. 

The uncertainty in the parameters can often be best de- 
scribed via a probabilistic uncertainty, e.g., 

A= {E{a;} = (w-;i)(a;- w)^} (11) 

where E {• • • } is with repsect to the underlying probabil- 
ity distribution of uj. If that is gaussian, then A = {w G 
N{Cj,R)}. 

Uncertainty also arises from noise in the control and other 
environmental sources, best represented by stochastic uncer- 
tainty, e.g., 

A = {^{uj{t)} = (I',E - Ld){Lu{t') - Cdf] = R{t, t')] 

(12) 

Although (9) is a convex set, as already mentioned, uncer- 
tainty modeling sets for SCP need not be convex, e.g., (11) 
and (12). Step 2 in (4) only requires that the uncertain param- 
eters be sampled from the uncertainty set. In our example to 
follow with uncertainty set ( 10) we use a simple uniform sam- 
pling. More sophisticated variants cycle though a sampling 
in the optimization step followed by validation on a different 
sampled set; bad parameters revealed in the validation step 
can be used in a new sampling for a repeat of the optimization 
step, e.g., [51]. 



N 



Ut=W exp 

fc=i 



T 



(8) 



C. Uncertainty modeling 

One general approach to modeling the uncertainty in the 
Hamiltonian parameters {uJz,lOx) is via the set describing a 
deterministic uncertainty. 



A = {\\W{u 



<1} 



(9) 



with parameter vector w = [oJx^'jJzY , nominal values a) = 
, oJzV , W a positive definite weighting matrix (here 2 x 2), 
and V typically 2 or 00. If = 00 and W is diagonal, then uj^ 
and are not correlated, in which case (9) is equivalent to. 



A _ / mill max min ^ 



<u,,<uj^-^} (10) 



If W is not diagonal, then uj^ and are correlated, possibly 
arising, respectively, from an approximation of a joint Gaus- 
sian or uniform distribution, with W, typically, the inverse of 



D. Robust optimization 

Gathering the previous pieces together we can form the fol- 
lowing instance of (3) for finding a robust control for a single 
qubit system by solving for the control magnitudes 9 e 
from, 

maximize min J^{0,uj3;,<^z) 

subject to Ut from (8) 

9 E Q from a combination of sets in Table I 
{uj^,uJa,) e A from (9), (1 1), or (12) 

(13) 

Since is a convex set and samples are taken from A to form 
the gradients (and possibly the Hessians), then Step 2 of (4) 
will be a convex optimization. 



IV. ROBUST ONE-QUBIT UNITARY OPERATIONS 

We use our SCP routine to construct control fields corre- 
sponding to the following set of one-qubit unitary operations 
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that are robust to variations in and Wz'- 



1 




1 1 




1 


1 


' ^/2 


1 -1 




e"/4 



(14) 



Here, Vi, Vh, Vp, correspond to identity, Hadamard, and 
phase gates, respectively. Note that Vh and Vp represent a 
universal gate set for one-qubit operations. The uncertainty 
set considered for all optimizations in this section is 

A - K e [0.99, 1.01], Lo, e [1.8, 2.2]} , (15) 

corresponding to 1% control amplitude uncertainty and 10% 
drift magnitude uncertainty. In each case, SCP is used to solve 
for 6* e from: 

maximize min T{d,ujx,uJz) 

= min |Tr(l/^C/T)/2|^ 

subject to Ut obtained from (8) 

9 e (unconstrained) 

(uJzjOJx) G A from (15) 

This is solved for all combinations of Vq, where a E 
{1, H,P}, and T e {1,2,4}, with select elements of iV e 
{5, 10, 20, 80}. For all of the results discussed in this section, 
the initial SCP controls c^"^ were optimized for unitary targets 
of the nominal Hamiltonian system 

H{t) = c{t)uj^X + ujzZ (17) 

(where uj^ is the average value of ijjx), such that 
J- [c^^^QxT^j^z] < 0.999, i.e., c^'^' is at best a nearly-optimal 
control field for the nominal system. Figure 1 presents the re- 
sulting control fields and corresponding fidelities. Properties 
of the robust controls are presented in table 11. There are eight 
configurations of N and T for each unitary target, represent- 
ing different combinations of N and T. With this one-qubit 
system and the uncertainty set A (15), controls with worst- 
case fidelities J^^c > 0.9999 are obtained for > 10 and 
T > 2 for all target operations. These results demonstrate 
that robust, high-fidelity control is possible with a relatively 
small number of control variables N, provided that the final 
time T is chosen properly. 

Interestingly, the worst-case fidelity for the identity gate ac- 
tually decreases as the number of time intervals increases from 
iV = 10 to iV = 20 for T = 2. This result suggests that the 
objective function in eq. (16) may possess local sub-optimal 
critical points, given that the set of controls with = 10 
and r = 2 is a proper subset of controls with = 20 and 



Target operation: Identity 



TV 


T 


logl()(l - -^wc) 


loglo(l --T^av) 


$[c] 


max[c] 


5 


1 


-3.13 


-3.82 


103.63 


16.21 


5 


2 


-2.35 


-3.16 


37.24 


7.25 


10 


1 


-3.28 


-4.20 


58.45 


11.88 


10 


2 


-5.23 


-5.79 


51.00 


6.59 


20 


1 


-3.31 


-4.24 


53.66 


13.47 


20 


2 


-4.35 


-4.98 


25.34 


6.03 


10 


4 


-4.62 


-5.66 


28.96 


4.22 


80 


4 


-5.08 


-5.60 


31.74 


6.00 



Target operation: Hadamard 



N 


T 


logio(l - -^wc) 




m 


max[c] 


5 


1 


-2.20 


-3.08 


32.27 


7.59 


5 


2 


-3.02 


-3.74 


16.35 


3.77 


10 


1 


-2.17 


-3.05 


36.93 


8.02 


10 


2 


-4.33 


-4.80 


30.33 


8.95 


20 


1 


-2.17 


-3.06 


34.38 


8.64 


20 


2 


-4.34 


-4.86 


30.07 


9.17 


10 


4 


-4.06 


-4.63 


16.60 


2.86 


80 


4 


-4.69 


-5.12 


25.61 


5.48 



Target operation: Phase 



N 


T 


logio(l - -^wc) 


loglo(l --^av) 




max[c] 


5 


1 


-2.77 


-3.51 


41.98 


8.39 


5 


1 


-3.71 


-4.19 


29.99 


6.70 


10 


1 


-2.96 


-3.55 


116.17 


28.62 


10 


2 


-4.34 


-4.88 


25.40 


6.80 


20 


1 


-3.02 


-3.61 


136.06 


33.88 


20 


2 


-4.30 


-4.77 


23.39 


7.12 


10 


4 


-5.57 


-6.02 


46.82 


5.87 


80 


4 


-6.00 


-6.34 


33.51 


5.91 



TABLE II: Properties of controls c{t, 6) for robust one-qubit oper- 
ations. Here, for the piece-wise-constant robust control from (13) 
<I>[c] = J^c^{t,6)dt = (T/N) \\e\\l, corresponding to control field 
fluence, and max[c] — \\0\\^ for maximum control magnitude. 



T ~ 2. For T = 1, the worst-case fidelity for the Hadamard 
gate also decreases as the number of time intervals increases 
from n = 5 to TV = 20. 



5 



(a) Control & fidelities for identity gate 





(c) Control & fidelities for phase gate 



N.5.T.1 N.5,T.2 N-IO.T-I N.10,T.2 N - 5, T - 1 N - 5, T - 2 N.10,T.1 N.10,T.2 




FIG. 1: Optimized control fields and corresponding fidelities for robustness to uncertainty {{uJx € [0.99, 1.01], ujx G [1.8,2.2]}) and for 
several different time intervals (A'^) and final times (T). (a) identity gate, (b) Hadamard gate, (c) phase gate. 
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V. GATE FIDELITY VERSUS CONTROL FLUENCE 

As a further illustration of the utility of SCP, we solve 
(13) with a constraint on control fluence and uncertainty in 
{ujx, ^z) in (5) for a target unitary of identity, V = /g, oper- 
ation time T ~ 2, and with the number of pulses N = 10. 
We solve this constrained optimization problem for a range 
of bounds on the fluence, denoted by 7, and for two uncer- 
tainty sets: A from (15) and a nominal case, Aq, effectively 
no uncertainty, 

Ao = K-2,c^, = 1} (18) 
Specifically, SCP is used to solve for 6 G R^" from: 
maximize min J^{9,Ldx,'^z) 

(lj^,lJj)gA 

= min |Tr(y^J7T)/2|^ 

subject to Ut obtained from (8) 

<P{9) = {T/N)\\0\\l < 72 for varying 7 
(w„cj,) e {Ao, A} from (18) and (15) 

For each uncertainty set {Aq, A} we start the SCP with the 
previously obtained robust control obtained from ( 1 3) for un- 
certainty set A (see the = 10, T = 2 plots in Fig. 1. The 
initial fluence constraint is set to 7 = 00, i.e., no control con- 
straint. We subsequently set 7 to 0.95 of the fluence of the re- 
sulting conti-ol ($(6*) = {T/N)\\e\\l with T/N = 0.2) found 
from the SCP routine and repeat the process, starting again 
with the previous initial control, reduced proportionally to the 
new fluence constraint, and so on. 

In Fig. 2, the soUd blue and dotted red lines show, re- 
spectively, the resulting tradeoffs between worst-case fidelity 
error log]^g(l — min(;^^ and achieved control flu- 

ence iT/N)\\e\\l for Ao (18) and A (15). The dashed blue 
plot shows the effect of the full uncertainty A on each of 
the respective optimized control fields obtained for the no- 
uncertainty set Aq as the fluence constraint is lowered. (Note 
that since we start with 7 — 00, i.e., no constraint, it follows 
that the right most points on each line correspond to uncon- 
strained controls.) The wiggly part of the solid blue plot for 
Aq is strictly due to the SCP stopping criterion: iterations halt 
if either log]^Q(l — < —16 or changes in log]^Q(l — T) are 
smaller than 10^^. It is important to note that it is not known 
if any of the points shown here are elements of a Pareto front 
for fidelity error versus fluence. 

Fig. 2 reveals some interesting, and initially surprising fea- 
tures. First, we see extreme sensitivity of the control opti- 
mized for the nominal parameter values Aq (solid blue) when 



T=2, N=10, =L 
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-14 




-16 



-a— e for with uncertainty 
- - - 9 for with uncertainty A 
■ 0- ■ ■ 9 for A with uncertainty A 



-a *i 



10 20 30 40 50 

Fluence, (T/N)||9||2 

FIG. 2: Worst-case fidelity error logj^Q(l — J-Wc) versus control flu- 
ence ^{e) = {T/N) \\e\\l from (19) for the identity gate Jg. For 
6 optimized for Ao (18), solid blue shows evaluation with Ao, and 
dashed blue shows evaluation with A (15). For d optimized for A, 
dotted red shows evaluation for A. 



evaluated over the full uncertainty range A (dashed blue). 
Specifically, at the right most points on the tradeoff curves, 
the latter control obtained from (19), denoted by 0a and the 
former, denoted by ^Aq^ are given by. 



7A = 



" -6.078 " 




" -6.075 " 


-6.557 




-6.554 


3.794 




3.798 


5.502 




5.505 


-1.756 




-1.760 


6.044 




6.040 


6.587 




6.583 


-3.753 




-3.757 


-5.507 




-5.511 


1.722 




1.725 



with the corresponding worst-case fidelities. 



min J'(6Iao) = 1 - lO-^^-^"" 

cjGAo 

min J"(6Iao) = 1 - 

min J'(6Ia) = 1 - lO^^ ^a 



1.000000 
0.999978 
0.999994 



Clearly the control optimized for A (dotted red) is only a bit 
better than the control optimized for Aq, and in addition, the 
control amplitudes are also extremely close. The main reason 
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for the closeness of the worst-case fidehties with respect to A, 
and the closeness of the control magnitudes, is that the start- 
up control, 6/^, gives very good performance. And a well- 
known phenomena with non-convex local optimizers is that 
solutions often tend to remain near the start-up when the start- 
up gives good performance [()8] which is certainly the case 
here. 

The extreme sensitivity to uncertainty is almost exclusively 
due to uncertainty in cj^,; variations in are not nearly as 
disruptive. This is due in part to the fact that in our exam- 
ple (5) ujx is a direct multiplicative uncertainty on the con- 
trol signal. Though small (1% variation), it is nonetheless 
in a sensitive spot. Gain variation in a classical open-loop 
control also exhibits this kind of sensitivity. Here, in the 
context of what is considered high performance for quantum 
computing, the sensitivity is extreme. To see the specific 
source we follow the procedure presented in First set 

uj^ — 1 + ujr^,\uJx\ < 0.01. Consequently, a variation in lu^ 
about Wj; = 1 is equivalent to a perturbation of 6* ^ 6 + uixO- 
This gives rise to the second order order Taylor series fidelity 
approximation. 

In our examples, the Hessian VgJ^(6', 1, cj^) is negative semi- 
definite, and in all cases the Hessian term dominates the gra- 
dient term. Although not shown, when evaluating the fi- 
dehty of any control over all uncertainty in ui^ and with 
\ujx I < 0.01, we get values of fidelity that coincide almost ex- 
actly with points on the tradeoff curve along the dashed blue 
plot. 

Another aspect revealed from the tradeoff curves is seen 
more readily in Fig. 3, which shows a few samples from the 
tradeoffs. Although visibly it is difficult to distinguish in 
Fig. 3(a) from 6'ao in Fig. 3(b), the fidelity surfaces are quite 
different. To some extent the opposite is true in Figs. 3(c-d) 
where the control magnitude differences are small but visible 
whereas the fidelity surfaces are almost the same, and clearly 
not so good for robustness. Table III summarizes some of the 
key features. In (a) the worst-case is slightly better than in(b) 
whereas in (b) the average is slightly better than in (a). How- 
ever, the somewhat central region in (b) where fidelity error 
is in the range of 10^^ is fairly large, suggesting that if the 
system (5) were designed differently we could have this much 
improved performance. So although in (b) the worst-case is 
worse than in (a) - although not so much worse - jointly con- 



(a) (b) 



0=51 N=10,T=2 *=50.99 N=10,T=2 
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FIG. 3: Robust control pulses and fidelity errors: (a) Start-up con- 
trol optimized for A (15) with $[c] < oo (no fluence constraint). 

(b) Control optimized for no uncertainty Aq (18) with $[c] < oo. 

(c) Control optimized for A (15) with <l>[c] < 33.83. (d) Control 
optimized for Aq (18) with $[c] < 32.57. 

sidering system design along with control can bring in unfore- 
seen possibly serendipitous designs. However, we also see 
from (c-d) that control constraints can quickly drive the de- 
sign out of the "sweet spot." Under such a fluence constraint 
would bring the design within reach to be one from Table II, 
say (for the identity gate) with T = 2 and = 20 fluence is 
25.43 which is within the constraints in Figs. 3(c-d). 



Control 
constraints 


Fluence 
achieved 


logio(l-.F) 


worst-case 


average 


best 


(a) A, $[c] < oo 


51.00 


-5.23 


-5.55 


-9.57 


(b) Ao,^c] < oo 


50.99 


-4.66 


-5.78 


-15.05 


(c) A,$[c] < 33.83 


33.83 


-2.65 


-3.20 


-6.94 


(d) Ao,<E>[c] < 32.57 


32.53 


-2.34 


-2.97 


-15.35 



TABLE III: Worst-case, average, and best-case fidelity errors for a 
few selected controls in Fig. 3 with fluence $[c] = (T/N) H^Hj- 



VI. SUMMARY 

Using SCP we have shown that a substantial fidelity ro- 
bustness is obtainable against uncertainties while simultane- 
ously using limited resources of control amplitude, bandwidth 
(number of pulses per time interval), and fluence. What is re- 
quired is a specific knowledge of the range and character of 
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the uncertainties, a process referred to in the control theory lit- 
erature as "uncertainty modeling." Although we have focused 
on a single example of a single qubit system, it is clear that 
even for this case there is a very rich variety of design trade- 
offs between control constraints (e.g., number of pulses, time 
interval, fluence, magnitude) and the range of system parame- 
ters from the underlying physics (e.g., the parameters {ujx,i^x) 
in (5)). Here we examined a 1 % variation due to lu^, which is a 
multiplicative gain perturbation to the control, and a 10% vari- 
ation in ujz, which to some extent is a coarse approximation to 
environmental uncertainty. The results presented give an indi- 
cation of what is possible for different values and ranges. For 
some implementations the physical parameters, and associ- 
ated ranges, can be "designed" by selecting the material used 
in the implementation, the circuit layout, etc.. Thus a com- 
bination of physical design and robust control design could 
lead to a "sweet spot" amongst the possibilities. One observa- 
tion which has surfaced in all our robust control examples is 
the sensitivity to uncertainty in the control implementation, in 
particular, the relatively small variation in ojx- This sensitiv- 
ity, though, should be viewed in the context of Fig. 2 where 
the fidelity error rises from essentially zero to a value in the 



range of 10~ . 

SCP is of course not the only approach to finding local so- 
lutions to non-convex problems. One advantage is the ease 
by which the design variable constraints and the uncertainty 
modeling can be directly incorporated in the local convex op- 
timization step of the algorithm. 

The array of results presented here hopefully herald what 
would be seen in more complex systems which involve mul- 
tiple qubits, the inclusion of controlled ancilla, the effect of 
coupling to a bath, and so on. In addition the results also be- 
gin to provide an insight into unanticipated control structures. 
Many of these potentialities are under consideration at present 
and will be forthcoming. 
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with modifications that take into account the specific underly- 
ing structure of the problem. 

Appendix B: Signal generation 



Appendix A: Convex optimization 

The convex optimization step in the SCP algorithm can be 
equivalently expressed as, 



maximize /o 

subject to fi + gje > fo, i = l,. 

e + 9eeje Otrust 



(Al) 



with f, = T{9, S^) e R and = WeJ'iO, 5.,) e R^. The 
optimization variables are now both ^ e R^ and /„. If both 
O and Otrust bound their respective elements in a "box" in 
R^, then (Al) is a linear program. 

The Hessian can be included by using its negative semidef- 
inite part, Ri = -—\^'^F{9^5i\- where [•]_ retains only the 
negative eigenvalues of the Hessian, all others set to zero. 
Then the worst-case fidelity constraint becomes, 

+ gJe ~ {1/2)9^ R,e > /o, i = 1, . . . , i 

Each of these inequlities is equivalent to a linear-matrix- 
inequality (LMI) in the variables (f?, /o) [ ]. The optimiza- 
tion step in SCP is now given by the semidefinite program. 



maximize /o 



subject to 



-k-gjo o^/V2 
e/V2 Rf 
ee, e e Gtrust 



> 0, i = 1, 



(A2) 



with Rf the pseudo-inverse of Ri. 

The optimization problems (Al) and (A2) are now in stan- 
dard forms suitable for use with existing software especially 
developed for these classes of convex optimization. In partic- 
ular, YALMIP [ ] and C VX [70, 7 1 ] are convex compilers 
compatible with MATLAB. Using these software tools makes 
it very easy to code the convex optimization problems almost 
exactly as expressed mathematically. These compilers call 
convex solvers such as SDPT-3 [ ] and SeDuMi [ ] which 
have been developed and in use for many years, and as a result 
are generally efficient and reliable. There are limits imposed 
by both memory and speed for a particular problem instance 
and computer platform. In these cases it is often necessary to 
use, if available, or devlop as necessary, specialized versions 



In general the control c{t, 9) is the output of a signal gen- 
eration device. As an example, consider the control generated 
from a system with time constant r and piece-wise-constant 
commands: 

C{t,e) - (1/T)(c(i,0)-C(i,0)), C(0)=0 

c{t, 9) = 6k, (fc - l)T/N < t < kT/N, fc = 1, . . . , iV 

(Bl) 

In this case. 



JV 



c{t, 9) = J2 Skit)9k = s{tf9, 0<t<T 



(B2) 



k=l 



This expression holds for any signal generation well rep- 
resented by known linear dynamics whose input is a se- 
quence of N control commands 9k,k = 1, . . . , at a uni- 
form sampling rate T/N. The linear dynamics are cap- 
tured in the shape function vector s{t) G R^, the elements 
of which in the above example (Bl) are given by Sfe(t) — 
1 _ g-{t-ik-i)T/N)/T fQj. _ i-^T/N < t < kT/N and 

Sfc(t) = (1 - e^T/N)/r^^^{t~kT/N)/r fo^ ^ > ^j./^^ j„ 

the ideal case with very fast dynamics (r — > oo) the con- 
trol is simply piece-wise-constant over N uniform time inter- 
vals of width T/N as given by (6) where the shape functions 
reduce to piece-wise-constant unit pulses of uniform width 
T/N. When the dynamics of the control generation device 
have an appreciable effect on these pulses, the product for- 
mula (8) for computing the final time unitary must be replaced 
by a continuous-time simulation. In our example this would 
be. 



iU = (c(i, 9)lu^X + uj,Z)U, Uq ^ I 



(B3) 



For any control in the form of (B2), the control constraint sets 
in Table 1 form convex sets in 9. For example, the fluence 
constraint <i>[c] = c^(t, 9)dt < 7^ is equivalent to ^{9) = 
9'^Qe < 72 with Q = /(f s{t)sit)'^dt. 

For laser control pulse shaping via a liquid crystal array, 
the controls are typically both phase and amplitude at each 
Fourier frequency component of the laser pulse after passing 
through an optical grating. Ignoring any dynamics in the crys- 
tal, the control signal is of the form c{t) — X^ili '^i sm{ujit + 
(f)i), with known frequencies oji and controls {ai,(j)i). The 
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control can be equivalently expressed in a form similar to (B2) 
as, 

K 

c{t,e)=Y,s,{tf6, (B4) 

i=l 

with multiple shape functions Si{t)'^ = [sm{ijjit) cos(a;ii] 
and Oj = ai[cos(j)i sin^^]. Since WOiW^ = ai, the magnitude 
constraint < < flmax brings (B4) again into a convex 



set in 9. The phases and amplitudes are easily obtained from 
any 9 so constrained found from SCR If, however, all the am- 
plitudes are the same, i.e., ai = gq, and only the phases are 
controlled, the constraint \\9i\\2 < ag is a convex relaxation 
[68] of the actual (non-convex) constraint ||6'i||2 = ao, and 
hence, SCP will return a local solution to the relaxed prob- 
lem. Some relaxations can be proven to be optimal; that is not 
known here, thus a post-optimization analysis is required. 
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